Astronomy <fc Astrophysics manuscript no. 3488 


February 5, 2008 


(DOI: will be inserted by hand later) 





A necessary extension of the surface flux transport model 

I. Baumann*, D. Schmitt, and M. Schiissler 



in 
o 
o 

(N 
-i— > 
O 

O 



> 

(N 

m 
O 

in 
o 

Oh 
i 

o 

o3 



13 



Max-Planck-Institut fiir Sonnensystemforschung, 37191 Katlenburg-Lindau, Germany 
February 5, 2008 

Abstract. Customary two-dimensional flux transport models for the evolution of the magnetic held at the solar 
surface do not account for the radial structure and the volume diffusion of the magnetic field. When considering 
the long-term evolution of magnetic flux, this omission can lead to an unrealistic long-term memory of the system 
and to the suppression of polar field reversals. In order to avoid such effects, we propose an extension of the 
flux transport model by a linear decay term derived consistently on the basis of the eigenmodes of the diffusion 
operator in a spherical shell. A decay rate for each eigenmode of the system is determined and applied to the 
corresponding surface part of the mode evolved in the flux transport model. The value of the volume diffusivity 
associated with this decay term can be estimated to be in the range 50-100 km 2 s -1 by considering the reversals 
of the polar fields in comparison of flux transport simulations with observations. We show that the decay term 
prohibits a secular drift of the polar field in the case of cycles of varying strength, like those exhibited by the 
historical sunspot record. 
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1. Introduction 



Sun: photosphere 



Flux transport models describe the evolution of the flux 
distribution at the solar surface as a result of the emer- 
gence of bipolar magnetic regions and the transport of 
the corresponding radial magnetic flux by the horizon- 
tal flows due to convecti on, differential rotation and 
merid i onal circulat i on (e.g . . lLeightonlll964t IbeVore et alJ 
19841 IWang et all [l989bt Ivan Ballegooiien et alJ 1998 : 
Schriiverll200li iMackav et alJ 120021: iDurrant et all 12004; 



jciuiiv^L]jzuu^| jMxicK.av eu au izuuzi mmianu eu ai.ii^uu'j 
Baumann et al.ll2004l). When app l ying t he m odel to mul 



tiple s olar cycles. ISchriiver et al.1 (|2002l) and lWang et alJ 

(2002) noticed that the secular evolution resulting from 
the simulations is not in agreement with the observations: 
activity cycles of varying strength lead to a drift of the 
polar fields and even to the disappearance of polar rever- 
sals. 

The conceptual deficiency in the conventional flux 
transport models that leads to this disagreement with the 
observations arises from ignoring the vectorial nature and 
the radial structure of the magnetic field and, particularly, 
from omitting the part of the diffusio n operator depend- 
ing on the radial coordinat e (see also I Wilson et al.l ll990: 
iDikpati fc Choudhurilll99l . This leads to an unwanted 
long-term memory of the surface flux since the decay time 
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of a dipolar surface field increases strongly in the presence 
of a poleward meridional flow. 

In order to avoid such unrealistic effects, IWang et alJ 
(2002]) made the assumption that the speed of the 
meridional flow is highe r du ring more active cy- 
cles, while ISchriiver et al.l (|2002) introduced an ad- 
hoc local decay term, calling upon the 'sea-serpent' 
flux emergence process combined with coron a l fiel d 
line reconnec tion as proposed by [Spruit et all l(l987lh 
Alternatively. IChoudhuri fc Dikpatil 1)19991 and earlier ref- 
erences therein) have considered a consistent MHD model 
for the evolution of an axisymmetric field in the meridional 
plane of the Sun. Unfortunately, our understanding of the 
solar dynamo and the available computational resources 
are still insufficient to carry out such simulations in the 
three-dimensional case. We therefore suggest to extend the 
flux transpor t model by a modal v ersion of the decay term 
proposed bv lSchriiver et alJ l)2002f) . We derive this term by 
considering the decay of the eigenmodes of the (volume) 
diffusion operator in a spherical shell. This leads to a con- 
sistent decay rate for each eigenmode, which can simply 
be incorporated in a code based upon expansion of the 
magnetic field into surface harmonics. The volume diffu- 
sion coefficient remains as a free parameter, which can be 
estimated by comparison of simulation results with obser- 
vational constraints (like polar field reversals) 

The outline of this paper is as follows. In Sect. [21 we 
determine the decay modes in a spherical shell for ap- 
propriate boundary conditions. The results are used in 
Sect.|3to define the decay term extending the surface flux 
transport model. The value of the volume diffusivity (free 



2 



I. Baumann et al.: A necessary extension of the surface flux transport model 



parameter in the deacy term) is estimated by consider- 
ing the evolution of the polar field over activity cycles of 
varying strength in Sect. 0] and by comparing the times 
of polar reversals with observational data in Sect. |S] We 
summarize our results in Sect.|SJ 



2. Decay modes of a poloidal field in a spherical 
shell 

In the abscence of systematic flows, the time evolution of 
a magnetic field is described by the diffusion equation, 



<9B 
~dt 



-77 V x (V x B) . 



(1) 



where 77 is the (constant) magnetic diffusivity. In spherical 
geometry, the magnetic field vector can be written as the 
sum of a poloidal and a toroidal part, the latter having 
no radial component. Since the flux transport model as- 
sumes a purely radial magnetic field at the solar surface, 
it is sufficient to consider the decay of a poloidal magnetic 
field in a spherical shell (representing the solar convection 
zone), setting the non-radial components to zero at the 
outer boundary (the solar surface). The evolution of this 
field is described by Eq. (Tj^), where B now represents a 
poloidal field. We introduce a spherical polar coordinate 
system with coordinates (r, 9, 4>), whose origin is located in 
the center of a sphere of radius Rq. The poloidal magnetic 
field can be represented by a scalar function, S(r,9,(j),t), 
as 



B = - V x (r x VS) = -r A S - 



d 
or 



(2) 



where A is the Laplace oper ator in spherical coordi- 
nates and r is the radius vector llBullard fc aellmanlll954t 
iKrause fc Radlerl 1 9801 . Inserting this field representation 
into Eq. (Tj^) we obtain 



V x r x V 



dS 
~dt 



= 7?VxVxVx(rx VS) 

= i)VxVx(rx AS) 

= -77 V x (r x VAS) , (3) 



where we have used V x r = 0. From Eq. © it follows 
that 



r x V 



77AS- 



dS_ 
~dt 



= 0. 



(4) 



The vector potential, — r x VS, is invariant under gauge 
transformations, so that S can be choosen such that the 
normalization condition 



1 Z7T 

J J S d(cos 0) 
-1 



d<j> = 



(5) 



is fulfilled for each value of r. Denoting by fl the angular 
part of the spherical Laplace operator, 



AS = 



1 



d_ 

Or 

1 d 
^888 



0_ 

Or 



,dS 



sin ( 



} dS 
~d0 



1 d 2 s 



sin t 



,dS 



we have 

1 2ir 



J J flSd{cos( 



= . 



(6) 



(7) 



With the normalization condition, Eq. (JSJ, we then find 

1 2tt 

AS d(cos 6) d<\> = . (8) 



-1 



Using this result, we obtain from Eq. (0} the scalar diffu- 
sion equation 



(9) 



We search for solutions of this equation within a spher- 
ical shell, ?*b < r < Rq, which are uniquely determined 
by specifying boundary conditions at the inner and outer 
boundaries. In accordance with the flux transport model, 
we assume the magnetic field to be radial at the surface. 
This leads to the boundary condition 



d(rS) 
dr 







at 



R 







(10) 



The bottom of the convection zone borders on the ra- 
diative core, which we represent by a field-free ideal con- 
ductor. Consequently, we require as boundary condition 
at r — rb that the radial component of the magnetic field 
vanishes, which is equivalent to the condition 



5 = 



at r = Tb 



(11) 



In the numerical calculations below we take rb = 0.7i?Q, 
which corresponds to the bottom of the solar convection 
zone. 

The general solution of Eq. © can be written as a de - 
compostion into orthogonal decay modes l)Elsasserlll94fj|) . 



00 oc 



n=0 1=1 m= 



where we have omitted the monopole term (I — 0). 
The functions Yi m are the spherical surface harmonics. 
Separation of variables leads to an exponential time de- 
pendence, viz. 



T ln (t) = exp(-?7 kf n t) , 



(13) 
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Table 1. Decay times r/„ = l/(r]kf n ) (in years) for a volume diffusion coefficient of 77 = f00km 2 s . The modes are 
characterized by I, the number of node circles on spherical surfaces, and n, the number of nodes in the radial direction. 





1 


2 


3 


4 


I 
5 


6 


7 


8 


9 





5.15 


4.42 


3.65 


2.97 


2.40 


1.96 


1.61 


1.34 


1.13 


1 


0.62 


0.60 


0.58 


0.56 


0.53 


0.50 


0.47 


0.44 


0.41 


n 2 


0.22 


0.22 


0.22 


0.22 


0.21 


0.21 


0.20 


0.20 


0.19 


3 


0.11 


0.11 


0.11 


0.11 


0.11 


0.11 


0.11 


0.11 


0.10 


4 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 



where 1 /i] kf n is the decay time of the mode characterized 
by the wave numbers I and n. For the spherical harmonics 
we have 



A Yu 



1(1 + 1) 



R% 



Y 



lm \ 



(14) 



so that we obtain the following differential equation for 
the functions Ri n (r): 



1 <l 2 Ri n 
dr 2 



2r 



dr 



"In 1 



1(1 + 1) 



Rln = . 



(15) 



Solutions are the spherical Bessel functions of the first and 
second kind, ji and yi, respectively, so that we have the 
general solution 



Rin(r) = a ln ji(k in r) + b ln yi(k in r) . 



(16) 



The linearity of the diffusion equation allows us to set 
the coefficients a;„ to unity without loss of generality. 
The lower boundary condition, Eq. (|ll|l . is then used 
to determine bi n = -ji(ki n r h )/yi(ki n r h ). Inserting this 
into Eq. I|16(l and using the upper boundary condition, 
Eq. (pi, leads to 



I 



ji(k ln R Q ) yi(k ln r h ) - ji(k ln r h ) yi(k hl R Q ) 
-ki n R® ji-i(k ln R Q ) yi(k hl r h ) 

-ji(ki„r h ) yi-i(ki n Ro) =0, 



(17) 



from which the eigenvalues fc;„ can be determined numer- 
ically. The index n is the number of nodes of the eigen- 
function Ri n (r) in r-b < r < Rq. Note that the above 
construction shows that the eigenvalues, fc;„, and, there- 
fore, the temporal decay rates of the eigenmodes do not 
depend on the azimuthal wave number, m, which describes 
the ^-dependence of the mode. 

Table^gives the decay times, n n = l/(r] kf n ), in years 
for various modes, assuming a volume diffusion coefficient 
of 77 = 100 km 2 s^ 1 . The dipole mode (I — 1) with n = 
has the longest decay time of about 5 years. While the 
decay times decrease only slowly for the higher multipoles 
(increasing I), the decay of the higher radial modes is much 
more rapid for the higher radial modes, being faster by a 
factor eight already for the mode (I = l,n = 1). ft is 
therefore justified to consider only the most slowly decay- 
ing modes with n = in the decay term for the surface 
transport model. 



3. Extension of the surface flux transport model 



We use the results obtained in Sect.|2to extend the surface 
transport model ijDeVore et al.lll98a IWang et al.lll989a|) 
by a decay term describing the volume diffusion of the 
poloidal magnetic field in the convection zone. This is cer- 
tainly only a rough description of the evolution of the sub- 
surface field, which can be improved once we have more 
quantitative information about the systematic flows in the 
deep convection zone. The working of the solar dynamo, 
on the other hand, is already represented in the model by 
the emergence of bipolar magnetic regions. 

The extended equation for the evolution of the surface 
flux is written as 



dB r 
~dt 



M0) 



dB r 

d<j> 



1 



d 



R Q sm6 86 



v(9)B r sin( 



Re 



1 d 



sin 6 86 



smt 



} 8Br 

86 



sin 



1 8 2 B r 



8(j) 2 



+ Q(e,4>,t)-D( v ) 



(18) 



where oj(6) is the angular velocity of the photospheric 
plasma, v(6) is the meridional How velocity on the solar 
surface, Q(6,<f>,t) is a source term describing the emer- 
gence of new magnetic flux, and rj^ is the turbulent mag- 
netic diffusivity associated with the nonstationary super- 
granular motions on the surface. We specify the decay 
term D(rf) on the basis of the decay modes in a spherical 
shell as determined in the previous section. To this end, 
we expand the instantaneous radial surface magnetic field 
into spherical harmonics, 



B r (R Q ,9,<p,t)=}^ 



00 m—-\-l 



Clm(t) Y im (0, (f>) 



(19) 



l— 1 m- 



The volume diffusion leads to exponential decay of each 
of these modes with the corresponding decay times, T\ n = 
1/(77 /c 2 n ), depending on the radial structure of the mag- 
netic field. Since the latter is unknown in the framework 
of the flux transport model, wc only consider the mode 
n = 0. This is justified by the fact that all higher modes 
(n > 1) decay much more rapidly (see Tabled, so that 
they do not affect the long-term, large-scale behaviour of 
the surface field that the flux transport models aim to 
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Fig. 1. Evolution of the north (solid) and south polar fields (dashed) in the case of a systematic variation of strength 
between even and odd cycles for four values of rj. Finite values of 77 lead to stable oscillations of the polar fields while 
the field reversals (zero crossings) occur earlier with respect to the preceding activity minima (indicated by the dotted 
vertical lines.) 



describe. We therefore write D(r]) as 

00 rn— +/ , x 

fl fa)=E E ^iwm). (20) 

l—l m— — l 

As we have seen in Sect. |2 the decay times do not de- 
pend on the azimuthal wave number, m, of the mode. The 
modal decay term given by Eq. i12()ll is particularly simple 
to implement in flux transport code s based upon an expan- 
sion i nto surface harmonics (e.g .. Ivan B allegooiien et alJ 
Il998l iMackav et al.ll20o3 iBaumann et al.ll2004j) . 

The decay term, D{rj), depends on the turbulent vol- 
ume diffusivity, r/, which will generally differ from the dif- 
fusivity 77h related to the flux transport by the horizontal 
surface motions of supergranulation. Since the velocities 
and length scales of the dominant convective motions in 
the deep convection zone, which are relevant for 77, are 
not well known, in the following sections we empirically 



estimate the value of r\ by comparing with observed prop- 
erties of the solar polar fields, namely the times of polar 
reversals and the sustained polar reversals in the case of 
solar cycles of varying strength. 

4. Cycles of varying strength 

One of the main problems with the original formulation 
of the flux transport model is the too long memory of the 
system. In the case of cycles with varying strength this can 
lead to a secular drift of the polar fields and a suppression 
of polar reversals. In the case of random fluctuations of the 
cycle amplitude, it results in a random walk of the polar 
fields superposed upon the cyclical variation. We illustrate 
this effects and their elimination through the decay term 
by way of a couple of examples. 

As a first illustration of the effect of the decay 
term introduced in the flux transport model, we con- 
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Fig. 2. Evolution of the north polar field (average above 75° latitude) for a sequence of synthetic cycles with the flux 
emergence rate taken proportional to the sunspot numbers since 1700 (full line) and 1750 (dotted line), respectively, 
for different values of the volume diffusivity, 77. The top panel (j]=0) shows the strong drift of the polar field resulting 
from the secular increase of solar activity during the last century (see also ISchriiverll200l|) . Finite values of rj reduce 
the secular effect and lead to more symmetric variations, as observed; the volume diffusion term also largely removes 
the dependence of the results on the initial conditions (difference between full and dotted lines). 



sider simulations of s ynthetic sets of activity cycles (see 
iBaumann et al.l [2004, for a detailed description of the 
model) that vary systematically in strength: the cycle am- 
plitudes (total amount of emergent flux) alternate by a 
factor of 2 between odd and even cycles. The results for 
the polar fields for several values of 77 are shown in Fig. ^ 
We define the polar field, -B po ie, as the averaged radial field 
strength over the polar cap poleward of 75° latitude. For 
?y = 0, i.e., in the absence of the added decay term, the po- 
lar fields show a strong secular drift (Fig. QJi) . This results 
from the fact that, during the weaker cycles, the amount 
of opposite-polarity flux reaching the poles is insufficient 
to cancel the existing polar field and to build up a field of 
opposite polarity. This results in a systematic drift of the 
polar fields. Including the decay term leads to a shorter 



memory of the system and a stabilization of the oscillation 
(with some fluctuations due to the random component in 
the prescription of the flux emergence). For growing val- 
ues of 77, the drift of the polar fields ceases earlier and the 
asymmetry as well as the amplitude of the oscillation de- 
creases (Fig. ^3-d). At the same time, the sign reversals 
of the polar fields occur earlier after the minima of flux 
eruption (activity minima, indicated by the dotted verti- 
cal lines) because the polar fields from the previous cycle 
have already been reduced by the effect of the decay term. 

We have also considered the secular variation of solar 
acti vity in the histo rical record of sunspot numbers (see 
also ISchriiveil l200l|) . Using the record of sunspot num- 
bers since 1700, we have simulated a series of solar cycles 
with the flux transport code. We have used the 'stan- 
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dard' parameters (butterfly diagram of emerging bipo- 
lar regions, tilt angles, polarity ru les, transport parame- 
ters, etc.) of lBaumann et all 120 041. except for taking the 
emergence rate of the bipolar regions proportional to the 
sunspot numbers (monthly sunspot numbers since 1750 
and monthly values interpolated from the yearly sunspot 
numbers before) , thereby also using cycle lengths in agree- 
ment with the observed record. The overall free scaling 
factor for the flux emergence has been fixed by requiring 
that the total (unsigned) surface flux during the last three 
solar maxima (cycles 21-23) for the c ase 77 = 100 km 2 s - 1 
matches the observed values given bv lAree et all l|2002|) . 
The evolution of the polar fields according to simulations 
with different values of 77 is shown in Fig. PJ1 In order to 
evaluate the effect of the arbitrary initial condition (zero 
surface field) we also show a run starting from 1750, rep- 
resented by the dotted line. In the case 77 = (no vol- 
ume diffusion), the long memory of the system leads to 
a drift of the polar fields in the 20 th century due to the 
secular increase of solar activity, so that the oscillations 
become very asymmetric, in striking contrast to the ob- 
served evolution of the polar fields. Finite values for 77 of 
the order 100 km 2 s _1 lead to more symmetric oscillations 
and a suppression of the unrealistic drift. The amplitude 
of the simulated polar field in the last cycles for a value 
of 77 = 100 km 2 s _1 is roughly consistent with the pub- 
lished observational data, which i ndicate amplitudes for 
the field strength of 10-20 G ^.g.-IWang fc Sheeleylll995l 
lArge et alJ2002HDikpati et al.l2004l:lDurrant et all2004j) 
Comparing the two runs with different starting times (full 
line and dotted line) , we find that the long memory of the 
system in the case 77 = leads to a significant difference 
between these two runs. For rj ^ 0, on the other hand, 
there is almost no dependence on the initial condition af- 
ter a few cycles. 

5. Polar field reversal times 

Having shown that the decay term prohibits drifts and an 
unrealistic long-term memory of the polar fields for values 
of 77 of the order of 100 km 2 s _1 , we now consider its effect 
on the calculated reversal times of the polar fields and 
compare quantitatively with observational results. 

5.1. Synthetic cycles 

In order to illustrate the effect of the decay term on the re- 
versal times, we consider synthetic cycles of equal strength 
and determine the reversal times of the polar field in de- 
pendence on the value of r\. In addition, we also consider 
the averaged (unsigned) field over the whole surface, -Btot, 
and the maximum polar field during a cycle. The abso- 
lute values of the field strength are arbitrary; here we are 
only concerned with the dependence on 77. Fig. [3] shows: 
i?tot for cycle maxima and minima, respectively (Fig.|3|i), 
the maximum polar field (Fig. 0t>), and the reversal time 
in years after the previous flux emergence (activity) min- 
imum (Fig. 0t), all as functions of 77. 
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Fig. 3. Synthetic solar cycles: dependence on 77 of the (a) 
averaged unsigned surface field (maxima and minima, re- 
spectively), (b) the maxima of the polar polar field, and 
(c) the polar reversal times after the previous activity min- 
imum (minimum of the cyclic flux emergence rate in the 
model). 

While the average surface field during cycle maxima 
varies only slightly with 77, there is a somewhat stronger 
decline of the values during cycle minima (Fig.[2la) , similar 
in proportion to the decline of the maximum polar field 
with increasing 77 (Fig.rjb). This results from the fact that, 
during activity minimum, both the polar and the total 
field are dominated by the dipole mode, which suffers from 
enhanced decay due to volume diffusion. A strong effect 
of the value of 77 on the polar reversal times is clearly 
visible in the lower panel. For increasing 77, the reversals 
occur earlier after cycle minimum, varying between 5.1 
years for 77 = and 2.8 years for 77 = 200km 2 s~ 1 . This 
results from the stronger decay of the polar field from the 
previous cycle, so that less opposite-polarity flux of the 
new cycle is required to reverse the field, leading to an 
earlier reversal time. 

Although these synthetic cycles are not intended for a 
quantitative comparison with actual solar data since they 
do not include a long-term modulation of the cycle ampli- 
tudes and also no variation of the cycle length, we may 
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Table 2. Epochs of reversals (average of north and south pole) of the simulated polar magnetic field poleward of 
±75° latitu de for three values of t he magnetic volume diffusivity, rj (in km 2 s _1 ), in comparison with the reversal times 
inferred bv lMakarov et ail l|2003h on the basis of polar crown filaments, T c , and from Ha synoptic charts, T\—i. The 
row at the bottom gives the average time intervals (in years) between the reversals and the preceding sunspot minima 
(Imirn last column). 
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nevertheless ch eck whether they are roughly consistent 
with the data. iMakarov et alJ l|2003h used two different 
methods to determine polar reversal times for the period 
1878-2001: (1) by the disappearance of the polar crown fil- 
aments, and (2) from estimating the magnetic dipole con- 
figuration on the basis of Ha synoptic charts. They find 
that reversals occur on average 5.8 years after the previous 
sunspot minimum for method (1) while they obtained a 
value of 3.3 years with method f2V IMakarov et all (2003) 
consider the first method to be more reliable. On the other 
hand, the magnetograph data compiled by lArge et all 
l|2002|) and iDikpati et all l)2004h for the the solar cycles 
21-23 indicate an average reversal time about 4.6 years 
after sun spot minimum, about the mean of the two values 
given bv IMakarov et all |2003) . Comparing this with the 
bottom panel of Fig. |21 we find that a value of ij in the 
range 50-100 km 2 s _1 leads to reversal times that are con- 
sistent with these results. A quantitatively more reliable 
comparison is carried out in the subsequent section. 

5.2. Simulation of solar cycles No. 13 - 23 on the 
basis of RGO/SOON sunspot data 

For a quantitative comparison with observations, we de- 
termine the magnetic flux input into the flux transport 
model on the basis of the sunspot group areas from the dig- 
itized version of the Royal Greenwich Observatory (RGO) 
photographic results, which are available for the time pe- 
riod 1874-1976. For the time after 1976, this series is con- 
tinued by data from the Solar Optical Observing Network 
(SOON) of the US Air Force. We have combined both 
datasets and transformed them into a sequence of emerg- 
ing active regions. Below we give a brief summary of the 
procedure; a more detailed description will be given else- 
where. 

We take each observed sunspot group to provide a 
bipolar magnetic region and determine its total unsigned 
magnetic flux from the observed sunspot area and an as- 



sociated facular (plag e) area according to the empirical re- 
lationship derived bv lChapman et all l)l997|) . We use the 
sunspot groups in the databases at the time of their maxi- 
mum area development, thus ensuring that every group is 
considered only once. Furthermore, we take only sunspot 
groups that reach their maximum area within ±45° of the 
central meridian. Owing to the large number of observa- 
tions, we consider this 90° window to be representative 
for the whole sun. In order to roughly account for the flux 
emergence on the full solar surface, we copy this 90° win- 
dow three times into the remaining 270° longitude range. 
Each of the sunspot groups determined in this way pro- 
vides a bipolar magnetic region whose orientation of the 
magnetic polarities is assigned according to Hales' polar- 
ity rules and whose tilt angle, a, is assumed to be given by 
Joy's law, i.e. a = 0.5 A, where A is the latitude. Further 
details a bout the treatment of b ipolar magnetic region are 
given in lBaumann et all ((2004) . The overall scaling factor 
in the magnetic flux of the emerging bipolar regions is 
chosen such that the flux transport simulations reproduce 
the evoluti o n of t he observed total surface field given by 
Arge et al. (2002) and the polar field strengths given by 
Dikpati et alJ l)2004j) . 

We have carried out a flux transport simulation based 
upon the set of input data covering the period 1874-2005. 
For the interval 1974-2005, the time-latitude diagrams of 
the longitudinally averaged magnetic field from the simu- 
lation can be compared with the corresponding diagram 
derived by D. Hathaway (NASA/Marshall Space Flight 
Center) from the NSO/Kitt Peak synoptic maps. Fig. 0] 
shows this comparison for three values of the volume diffu- 
sivity: r\ — 0, 100, 200km 2 s~ 1 . The evolution of the global 
field is best reproduced by the case i] = 100 km 2 s _1 , which 
yields the nearest agreement of the reversal times of the 
polar fields for the last three cycles. 

In order to compare with the result of IMakarov et all 
(2003), we have determined the epochs of the reversals 
of the average field poleward of ±75° latitude from the 
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Fig. 4. Comparison between observed magnetic field distribution from 1974 on (panel a, based upon NSO/Kitt Peak 
data, courtesy D. Hathaway) and flux transport simulations for different values of the volume diffusivity: r\ = (panel 
b), i) = 100km 2 s -1 (panel c), and r\ = 200km 2 s~ 1 (panel d). Shown are time-latitude diagrams of the longitudinally 
averaged surface field. The simulations reproduce the global development of the magnetic field, with poleward surges 
of following-polarity flux leading to reversal of the polar fields around sunspot maxima. However, in the case without 
the decay term (panel b), the reversals occur too late. For non- vanishing values of r], the reversal times are shifted 
forward in time; the best agreement with the observed reversal times is obtained for rj = 100km 2 s _1 (panel c). 
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flux transport simulations for the solar cycles 13-23. The 
reversal times (average of the north and south polar re- 
versals) are given in Table 01 together with the reversal 
times inferred by iMakarov et all l)2003[) from the evolu- 
tion of polar-crown filaments (T c ) and from Ha synoptic 
charts (Xi = i), respectively. The last row gives the average 
time interval (in years) between the polar reversal and the 
preceding solar minimum. We recover the trend towards 
earlier reversals for larger values of i] found in the pre- 
vious section for the synthetic cycles. The values cover 
about the ra nge defined by the resu lts from the two meth- 
ods used bv IMakarov et all ((2003) . However, even with 
rj = 50 km 2 s _1 we do not reach the value of 5.8 years that 
they determined from the disappearance of polar crown 
filaments. This could be due to the fact that we define the 
polar reversal from the average field in the caps poleward 
of ±75° latitude, while the polar crown filaments probably 
disappear later, when the last remnant of the old polar- 
ity in fact vanishes. A fair degree of avera ging is also in- 
herent in the magnetograph data shown bv lDiknati et alJ 
l|2004h . from which we roughly estimate an average rever- 
sal time of 4.6 years after sunspot minimum. Altogether, 
we conclude that a volume diffusivity in the range 50-100 
km 2 s _1 is probably adequate for flux transport models 
of the solar surface field. This is also consistent with the 
timings of the polar reversals with respect to the sunspot 
maxima: for r/ — 50 km 2 s _1 the polar field reverse on av- 
erage 1.3 years after solar maximum; for rj = 100km 2 s _1 
we find an average value of 0.3 years after maximum. 



6. Conclusion 

We have shown that a modified (modal) version of the 
de cay term first introdu ced into the flux transport model 
bv lSchriiver et al 1 l|2002l) can be derived consistently from 
the volume diffusion process that is (among other factors) 
neglected in the standard flux transport models. Including 
this term removes the unrealistic long memory of the sys- 
tem and thus prohibits a secular drift or a random walk of 
the polar fields in the case of activity cycles with variable 
amplitude. The value of the (turbulent) magnetic volume 
diffusivity, rj, can be estimated by considering the reversal 
times of the polar field relative to the preceding activ- 
ity minimum and by comparing with direct observations 
or values inferred from proxy data. We find that values 
of T) in the range 50-100 km 2 s _1 are consistent with the 
observational constraints. This is also the order of mag- 
nitude suggested by simple estimates based on mixing- 
length models of the convection zone. With this consistent 
extension and improvement of the model, flux transport 
simulations of the large-scale magnetic field on the solar 
surface over many activity cycles can be carried out. 
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